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Abstract. We set-up and solve the Cauchy problem for Schrodinger-type differential operators 
with generalized functions as coefficients, in particular, allowing for distributional coefficients in 
the principal part. Equations involving such kind of operators appeared in models of deep earth 
seismology. We prove existence and uniqueness of Colombeau generalized solutions and analyze 
the relations with classical and distributional solutions. Furthermore, we provide a construction 
of generalized initial values that may serve as square roots of arbitrary probability measures. 



1. Introduction 

Partial (and pseudo-) differential operators of Schrodinger-type with variable coefficients in the 
principal part arise in the form of so-called paraxial equations in models of wave propagation. 
These are based on narrow-angle symbol approximations of wave operators and have been applied 
in integrated optics, underwater acoustic tomography, reflection seismic imaging, time-reversal 
mirror experiments (cf. [2, 5, 41]), and recently in [12] to seismic wave propagation near the core- 
mantle boundary inside the Earth (at approximately 2800 km depth). Since paraxial equations 
are used to split the wave fields according to a prescribed principal direction of propagation they 
are also called one-way wave equations. The leading-order approximation leads to model equations 
of Schrodinger-type, where the material properties are still reflected by the regularity structure 
of the coefficients in the principal part. Under strong smoothness conditions on the wave speed 
function (i.e., the coefficients in the original wave operators) well-posedness of the one-way wave 
Cauchy problems has been discussed in [20, 39]. In [12] this smoothness assumption has been 
considerably relaxed by allowing the coefficients to be of Holder- or Sobolev-type regularity below 
log-Lipschitz continuity. Recall that in general existence of distributional solutions to the second- 
order wave equation may fail below log-Lipschitz regularity of the coefficients (cf. [10]). Beyond 
such coefficient regularity barrier unique solvability of the Cauchy problem holds in the sense of 
Colombeau generalized functions (cf. [17]), even in a covariant setting. 

The relevance of coefficients with low Holder or Sobolev regularity has been shown in a variety of 
geophysical applications, e.g. in the study of phase transitions in Earth's lowermost mantle ([40]) 
or in exploration geophysics ([13, 21, 22, 24, 32, 36, 42]). The model analyzed in [12] describes 
the paraxial approximation to seismic wave propagation near the core-mantle boundary by the 
following Schrodinger-type equation with depth z as evolution variable, the 2-dimensional lateral 
x-variable, and a pseudodifferential time-frequency dependence 

(1) dzu- i{dxi{cidxiu) + dx2{c2dx^u)) =0, 
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where Ci{z, x, r) and C2{z, x, r) are strictly positive symbols of order — 1 in r, continuously differen- 
tiable with respect to z, but a non-Lipschitz dependence regarding the .x-regularity, e.g. of Sobolev 
type H'^'^^ or Holder continuous with exponent r, where r < 1. Roughly speaking, the main result 
of [12] is the unique solvability of the corresponding Cauchy problem for the wave component u 
in C([0, oo[, c5^'(IR; _ff^(M^)) (i.e., continuous dependence on z, arbitrary ,5^' quality with respect 
to time or frequency, but lateral x-regularity of class if^), given initial data in L^(R, H'^iM?)) (at 
2 = 0) and a right-hand side in e^([0, oo[, L^(]R, L^(]R^))). The emphasis in this result is on the 
(Holder or) Sobolev regularity of the solution (with respect to x) in relation to the initial data 
regularity under lowest possible regularity assumptions on the coefficient (cf. the brief discussion 
on inverse analysis of medium regularity at the end of [12]). 

In the current paper we strive for a widening of the possible range of applications and thereby 
also simply for a more unified general mathematical set-up. Thus, we allow for discontinuous or 
distributional coefficients, initial data, and right-hand sides. In fact, a natural method is then to 
place the whole Cauchy problem into the context of nonlinear theories of generalized functions, 
in particular the differential algebras of Colombeau generalized functions (cf. [7, 8, 9, 19, 29]). In 
view of models from geophysics this amounts to including discontinuous material properties (e.g. 
at fault zones and geological boundaries or cracks) and Dirac-type data such as strong impulsive 
sources (e.g. explosions or earth quakes). A quantum field theoretic application of Colombeau 
generalized functions in terms of regularizations of Wightman distributions has been given in 
[18]. In the simpler context of quantum mechanics one will of course be interested in allowing 
for a singular zero order term in the operator to represents a potential and non-smooth data 
corresponding to an initial probability distribution. The square of the modulus of a solution 
to the standard Schrodinger equation is usually interpreted as (time evolution of a) probability 
density. Why should one not turn this around in the sense of allowing now for generalized initial 
data which represent a "square root of a given arbitrary probability measure"?^ We will show 
below how to construct a Colombeau generalized function whose square is associated with a given 
probability measure in the sense of distributional shadows. The issue of squares of distribution 
theoretic objects or generalized functions which contain or model measures (in particular, Dirac- 
type terms) also arises in general relativity theory, e.g. in the pseudo-Riemannian metric associated 
with impulsive pp-waves (cf. [19, Section 5.3] and [23, 34, 35]) and in the stress energy tensor 
corresponding to ultrarelativistic Reissner-Nordstr0m fields (cf. [33]). A different regularization 
approach for powers of Dirac measures as initial value appears in [27] in the context of (semilinear) 
heat equations with singular potentials. 

In our analysis of the Cauchy problem we will return to denoting the evolution variable by t and 
furthermore incorporate also first-order terms in the differential operator. Note that unlike in the 
specific geophysical model situation of [12] we neglect additional pseudodifferential aspects in the 
symbol, since these played only the role of an external parameter upon a partial Fourier transform. 
In summary of the basic structure, let T > be arbitrary. We consider the Cauchy problem for a 
generalized function u on R." x [0, T] in the form 



where Cfc (fc = 1, . . . ,n), V , and / are generalized functions on M" x [0,T] and g is a general- 
ized function on M". We note that basics of a theory for abstract variational problems in the 
context of Colombeau-spaces have been presented in [16]. As indicated in [16, Section 8] the 
Lax-Milgram-type theorem established there provides solutions to Diriehlet-problems involving 
differential operators similar to the spatial derivatives appearing in the left-hand side of our above 
differential equation. However, this does not directly extend to our kind of evolution problem and 



n 



(2) 




(3) 



u |t=o = 9, 



My thanks go to Gebhard Griibl and Michael Oberguggenberger who brought up this viewpoint in joint 
discussions at the University of Innsbruck in April 2002. 
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furthermore involves less strict solution spaces than we are willing to accept here, namely Shi as 

compared to ■ Colombeau-generalized solutions to linear and nonlinear Schrodinger equations 
with constant coefficient principal part have been constructed previously in [4, 37, 38]. 

We may point out that our current paper has both focus and aim different from the enormous 
literature on the prominent questions concerning properties of Z/^-solutions and spectral theory for 
Schrodinger equations with standard principal part (i.e., with constant or smoooth coefficients) 
and singular potentials. Instead we investigate here the feasibility of extending the solution concept 
and the basic analysis to the case of generalized functions and, in particular, non-smooth or highly 
singular principal parts. The plan of the paper is as follows. In Section 2 we briefly review the 
regularization approach to generalized functions in the sense of Colombeau and show how square 
roots of probability measures can be implemented in this framework. Section 3 presents the main 
result of unique existence of generalized solutions to our Schrodinger-type Cauchy problem (2- 
3). Furthermore, we discuss the relation of Colombeau generalized solutions with classical and 
distributional solution concepts. 



2. Regularization of coefficients and data 

The basic idea of regularization methods is to replace non-smooth data by approximating nets 

of smooth functions, e.g. instead of g consider C°° 9 g {e ^ 0). More generally, we may 

replace ^ by a net {ge)o<e<i of functions, convergent or not, but with moderate asymptotics 
with respect to e and identify regularizing nets whose differences compared to the moderateness 
scale arc negligible. For a modern introduction to Colombeau algebras we refer to [19]. Here we 
will also make use of constructions and notations from [15]. 

Construction of generalized functions based on a locally convex topological vector space E: 

Let E he a. locally convex topological vector space whose topology is given by the family of 
seminorms {pj}jej- The elements of 

Me ■■= {{ue)e e E^°'^^ -.yj & J3N eN pj{ue) = 0{e-^) ase 0} 

and 

:= {{u,), e : Vj e J Vg e N pj{u,) = 0{e^) ase ^ 0}, 

are called E-moderate and E-negligible, respectively. With operations defined componentwise, e.g. 
(Ue) + {vs) := («£ + Vg) etc., 'Ne becomes a vector subspace of M^;. We define the generalized 
functions based on E as the factor space Se '■= Me/^e- If is a differential algebra then T^e is 
an ideal in M^; and 3e becomes a differential algebra too. 

By particular choices of E we reproduce several standard Colombeau algebras of generalized 
functions. For example, E = C with the absolute value as norm yields the generalized complex 
numbers Sb = C and for any 17 C M'' open, E = with the topology of compact uniform 

convergence of all derivatives provides the so-called special Colombeau algebra = S(fi)- Recall 
that J7 I— > S(f7) is a fine sheaf, thus, in particular, the restriction u\b of u <E S(f7) to an arbitrary 
open subset i? C is well-defined and yields u\b G S{B). Moreover, we may embed D'{^) into 
S(n) by appropriate localization and convolution regularization. 

In case E C D'(il) certain generalized functions can be projected into the space of distributions 

by taking weak limits: we say that u G Sb is associated with w G D'{^1)^ if ti^ ^ w in D'{VL) as 
£ — > holds for any (hence every) representative (wg) of u. This fact is also denoted by w « w. 

In the current paper we will consider open strips of the form Vlx = K"x ]0, T[C (with T > 

arbitrary) and employ the spaces E = H°°{nT) = {h e e°°{Tl^) : d°'h e L'^inr) Va e N"+i} 
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with the family of (semi-)norms 

\ce\<k 

as well as E = W°°^°°{Qt) = {h G e°°{Th) : d"h e L°°(Ot) Va e N"+i} with the family of 
(semi-)norms 

\\h\\^,,^ = m^\\d'^h\\^^ (fceN). 

|c«|<fe 

(Note that clearly satisfies the strong local Lipschitz property [1, Chapter IV, 4.6, p. 66] and 
therefore the Sobolev embedding theorem [1, Chapter V, Theorem 5.4, Part II, p. 98] implies that 
every element of iJ°°(f2T) and W°°'°°{nT) belongs to e°°{n^).) 

To avoid overloaded subscripts we shall make use of the following (abuses of) notation in the 
sequel 

3^2 (M" X [0,r]) Sn^in^) and S^^ (M" x [0,r]) 9w^...^n^y 

For example and in explicit terms we will represent a generalized initial value g € Sl^ (R" x [0, T]) 
by a net (gg) with the moderateness property 

Vfc3m: ll^^ll^, = 0(e-™) (e ^ 0) 

and similarly for the right-hand side and the coefficients. Asymptotically negligible errors in the 
initial data, e.g. by using a representative (5^) instead, are expressed by estimates of the form 

VfcVp: \\g, - g,\\jj, = 0{sP) {e ^ 0). 

Similar constructions and notations will be used in ease oi E = _ff°°(M") and E = 14^°°'°° (R"). 
Note that by Young's inequality ([14, Proposition 8.9. (a)]) any standard convolution regularization 
with a scaled mollifier of Schwartz class provides embeddings and Sl^ (1 < P < 

00). 

As an example of a detailed regularization model we construct Colombeau generalized positive 
square roots of arbitrary probability measures, which can serve as initial values in the Cauchy 
problem analyzed below. 

Proposition 2.1. Let jj, be a (Borel) probability measure onM". Choose p e Li(R")nVF°^^°^(M") 
to be positive with J p = 1 and satisfying p{x) > \x\~'^° when \x\ > 1 with some mo > n. Set 
pg(x) = -prp{^) o.T^d hg := p,* Ps, then we have 

(1) hg is positive and setting (f)^ := \/h^ the net {^s)ee]o,i] represents an element (j) e S(K"') 
such that (j)^ « p; 

(2) there exists a generalized function g G Si2(M") such that g^ ^ p and the class of {ge\Q.)ee]o,i\ 
is equal to (I>\q in 9{0,), or by slight abuse of notation g\Q = (p\a, for every bounded open 
subset new. 



Proof. (1) The Borel measure p is regular since it is finite and E" is locally compact and second 
countable ([6, Proposition 7.2.3] or also [14, Theorem 7.8]). Thus //(R") = 1 implies that we 
can find a compact subset A C R" such that p{A) > 1/2. A variant of Young's inequality for 
measures (cf. [14, Proposition 8.49]) applied to d^h^ = p * d"ps directly implies that the net 
(he) is C°°(R"')-moderate, even with global i°°-norms, since \\p* 9"pe||^oo < ll/^llvar ll^"PellL°° — 
II^^pIIloo^"""'"'- (Here || \\^^^ denotes the total variation norm on the space of finite Radon 
measures on R", which gives llA^llvar — p{^^) = 1 in case of the positive probability measure p). 

Furthermore, for any x G M" and e > we have 

he{x) = p * ps{x) = / pe{x - y)dp{y) > I pe{x - y)dp{y) > p{A) ■ min ps{z) > Q 
J J ze{x}-A 

R" A 
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and hence that cj)^ = y/h^ € e°°{W) {0 < e < 1). Moreover, if if C M" is compact then we have^ 
with r{K) := max{|a;| : x G K — A} > the estimate 

(*) inf he{x) = inf /x * Pe{x) > ■ inf min Pe{z) > - ■ min ^ ^ 



> — min p{z/e) > 7-177 — = ^ ,,,, — (0 < £ < 1/rm). 

If a S N" is arbitrary then 9"(^£ is a linear combination of terms of the form d^'^he ■ ■ ■ d^^hejlii'^ 
with appropriate /3i, . . . , l3k € N" and ? G N. Hcnco C^(R")-moderateness of (he) together with 
the lower bounds for infxeK he{x) obtained above prove C°°(M")-moderateness of <j>e- 

Finally, wc have that by construction 0^ is represented hy hg = fj, * Pe and thus clearly converges 

to /i as e — > in the sense of distributions. 

(2) For j e N let Kj be the closed ball of radius 2^ around in M". Let xo S 2)(M"), < xo < 1 
with xo{x) = 1 when |a;| < 1 and Xo{x) = when \x\ > 2 and put Xji^) •= Xo(2~-'a;) {j > 1). 
Thus we have Xj{^) = 1 when \x\ < 2^, Xji^) = when \x\ > 2-'+^, and for every 7 e N" 

We define the net {ge)e&]o,i\ of smooth functions on K" by 

g, := Xj ^e, if 2-^-' < e < {j € N). 

For every e G ]0, 1] the function is smooth and has compact support, hence belongs to 
We have to show that (ge) is an iI°°-moderate net. 

Let a G N" arbitrary and 2~^~^ < e < 2~^. By the Leibniz rule we have that d^g^ is a linear 
combination of terms of the form d^cf)^ ■ d^^^Xj iP ^ ck)- Therefore wc have by the triangle 
inequality that ||9".g£||^2 is bounded above by a linear combination of terms of the form 

\\dP^,-d-^Xjh^ < Wd'-^X.h^WdUeh^^K,^,) 

and it remains to prove e-moderate bounds for 'i'e^j^'yo i^x ■ ^i) when /3 < a (note that the coupling 
of j with e has to be taken into account here). As already noted in (1) the latter is in turn bounded 
above by a linear combination of terms of the form 

\d^^h,{x)\---\d^''h,{x)\ 

with suitable /3i,...,f3k G N" and I G N. Thanks to the estimate (*) in the proof of (1) we 
have l/|/i£(x)|'/^ < 2r{Kj+i)'^°e~''"^°~"-'> uniformly with respect to x G Kj+i, where r{Kj+i) = 
0(2-'+^) = 0(l/£). Thus it remains to consider the factors \d^''hs{x)\ in the above supremum. 

As has already been noted in the first paragraph of the proof of (1) we have global _L°°-estimates 
with moderate e-dcpendcncc. Thus in summary, we have shown the i/°"-modcratcness of (g^). 

Finally, let 17 C R" be bounded and open. Choose j sufficiently large so that Q C Kj. Then 
we have for < £ < by construction of g^ that = (j)^\n and hence equality of the 

corresponding classes in 9{^)- Moreover, the latter also implies that g^ « p,, since the support of 
any test function in D(R") is contained in some open bounded subset. □ 

Remark 2.2. (i) Property (*) established in course of the proof shows in fact that the class of 
(fte) in S(M") as well as cp are strictly positive generalized functions (cf. [26, Theorem 3.4]). In 
this sense, ^ provides a strictly positive square root of the positive measure p,. 



We may exclude the case A = K = {xq} without loss of generality 
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(ii) For specific choices of p in L^{W) n H°°{W) such that ^ e H°°{W) we could obtain that 

is also _ff°°-modcrate and directly defines a square root in 3^2 (R") without having to undergo 
the cut-off procedure in part (2) of Proposition 2.1 (which, on the other hand, cannot be avoided 
for general p G H°°). For example, putting p{x) = c(l + |a:;p)~^"+^^/^ with a suitable normalization 
constant c > provides such a mollificr. However, the above Proposition leaves considerably more 
flexibility in adapting the regularization to particular applications. 

3. Generalized function solutions and coherence properties 

We come now to the main existence and uniqueness result for generalized solutions to the Cauchy 
problem (2-3). We recall that a regularization of an arbitrary finite-order distribution which meets 
the log-type conditions on the coefficients Ck and V in the following statement is easily achieved 
by employing a re-scaled mollification process as described in [28]. 

Theorem 3.1. LetCk (k = 1, . . . ,n) as well as V be reaP generalized functions in Sloc(M"x [0,T]), 
/ in Si,2(R" X [0,T]), and g be in Si2(IR"). Assume that the following log-type conditions hold 
for some (hence all) representatives (cA:e)ee]o,i] of Ck (k = 1 . . . ,n) and (14) of V: 

WdtCkehoo = 0{log{^)) and WdtVehoo = 0{log{^)) (e ^ 0). 

In addition let the positivity conditions Cke{x,t) > cq for all {x,t) e M" x [0,T], e e]0, 1], k = 
1, . . . , n with some constant cq > 6e mei.^ 

Then the Cauchy problem 

n 

(4) dtu-i'^ a^fe (cfe^xfc u) - iVu = f, 

fe=i 

(5) u \t=o = g, 
has a unique solution u G Sl2(R" x [0,r]). 

Proof. In terms of representatives (4-5) means that we have a family of Cauchy problems parametrized 

by£e]0,l]: 

n 

(6) dtu, - i ^ 5,, {ckeda^^Ue) - iV,u, = f, on Ot = R"x ]0, T[, 

(7) ue |t=o = ge on K", 

where Cks {k = l,...,n) and Ve (are real-valued and) belong to W°°'°°{^It)^ cke satisfies the 
log-type and the positivity condition as stated above, g^ G iJ°°(IR"), and € H°°{Qt)- 

Our strategy is to solve the corresponding problem at fixed, but arbitrary, parameter value e 
and thereby produce a solution candidate (ug). The substance of the proof lies in the efforts to 
show that each belongs to if°°(M"x]0,T[) and in addition satisfies moderateness estimates in 
every derivative with respect to the L^-norm. We will obtain basic energy estimates by standard 
variational methods as discussed by Dautray-Lions in [11, Chapter XVIII, §7, Section 1], but we 
will need to perform an additional analysis of the dependence of the constants in these estimates 
on the various norms of the coefficient functions. 

Step 1 (basic estimaies and, regularity): We apply the set-up and constructions in [11, Chapter 
XVIII, §7, Section 1] to the Hilbert space triple H^{W) C L^{W) C i?-i(R") and the sesquilinear 
form 

n 

ae{t;v,^P) :=^(cfce(t)a,,¥>|a,,^)i2(K„) + (K(i)^|V) {ip,i; G H\R"j), 
fe=i 



'In the sense that they possess represcntating nets of real- valued functions. 

^For any other representative the conditions thus hold for sufficiently small £ with co/2 instead. 
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where we have used the short-hand notation v{t) for a function x v{x,t). The basic conditions 
[11, (7.1) and (7.2), p. 621] on the quadratic form are easily seen to be met: continuous 
differentiability of the map t i— >■ a£{t; (p, ip) follows from smoothness of c^e {k = 1, . . . ,n) and 14, 
the form as{t; ., .) is hermitian since each Cke and is real-valued, and the positivity conditions 
on the coefficients Cfee (fc = 1, . . . , n) in the statement of Theorem 3.1 yield the following coercivity 
estimate 

where Ag = Co+||T4||^oc (thus we have A = Ae and a = cq in [11, (7.2), ii) on p. 621] ). Furthermore, 
the hypotheses [11, (7.5) and (7.6), pp. 621-622] on the initial values hold by our assumptions on 
and fg. We claim that [11, Theorem 1, pp. 621-622] in combination with the remark on additional 
regularity [11, Remark 3, p. 625] implies that we obtain a unique solution 

(Si) ue e e([o,T],ifi(M"))nei([o,T],ff-i(M")) 

to the Cauchy problem (6-7) for every e G]0,1]. In fact, the basic theorem gives existence and 
uniqueness with G L^qq, T], ff ^(R")) and < G L'^{[Q,T],H-^{W')) only. However, as men- 
tioned in [11, Remark 3, p. 625] the existence proof given in Dautray-Lions directly shows that the 
solution as well as its t-dcrivativc are L°°-functions of t with values in or H^^, respectively. 
The same remark states that even continuity with respect to t holds (an earlier reference for this 
fact is [25, Remark 10.2, p. 302], which also indicates a proof based on convolution regulariza- 
tion). Finally, we note that thanks to e([0, T],H'^) C e([0, T], i?"^) an apphcation of [11, Chapter 
XVIII, §1, Section 2, Proposition 7, p. 477] shows that the solution belongs to ei([0, T], if-i). 

Moreover, a careful inspection of all constants appearing in the derivation of the a priori estimate 
[11, (7.16) and (7.17) on p. 624] shows that we may deduce the following basic energy estimate 
uniformly for all t G [0, T] 

(8) |k.(i)||Li <C2ee^- 1^115.11^1+ /(ll/.WIli^ + ll^tAMIlL-O'^^ 

where the specifics on the constants Ci^ and Cie according to the derivation in [11, pp. 623-624] 
are as follows: C2e = 0{T{co + HV^eH^oo)) = 0{e~P) {e 0) for some p G N, and 

Cie = o(- • ( max WdtCkehoo + WdtVeW^^)) {e ^ 0), 

VCq l<k<n / 

which by the log-type conditions on (cke) and (K) imply moderateness estimates for the first-order 
spatial derivatives in the form 

BNiGN: \\d^ML2mn^,0T])<'^ ^^^P \Kmm = Oie-""') {e ^ 0, I = 1, . . . ,n). 

^ " t6[0,T] 

As a preparation for the procedure in Step 2 we claim that additional regularity properties for 
hold in the form 

(*) u, G L°°{[0,T],H'^{W))nW^'°°{[0,T],H\R'')) 

(at fixed e without precise asymptotic estimates of the norms) . These can be obtained from the 
concept of mild solutions and evolution systems on L^(M") (cf. [30, Chapter 5]) for the self- adjoint 
familiy of operators {As{t))teio,T]> where Ag{t)v{x) := ^l^i dxk{cks{x,t)da:^v){x) {v G //"^(R")) 
with common (i.e., t-independent) domain iJ^(K"): In fact, Equation (4) has the form 

dtUe - iAe{t)Ue = fs + iVeUe =: Fe, 

where we already know from (5*1) that G e([0, T], ^^(K")) n &{[0,T], H-^iR"-)). Let the 
strongly continuous evolution system corresponding to {Ae{t))t^[o^Tl be denoted by (Ue{t, s))o<s<t<T, 
then we necessarily have by Duhamel's formula 

t 

u,{t) = U,{t,0)g, + j U,{t,T)F,{T)dT G ei,([0,r],L2(R»)), 



where C^([0,T],L^(IR")) means weakly continuously differentiable as t-dependent distribution on 
M" with values in L^{W) for every t. Writing the differential equation now in the form 

we may employ elliptic regularity (spatial, with respect to x) for the second-order operator Ag{t) 
and deduce G e^([0, T], iJ2(M")) C L°°([0, T], iJ2(R")). Moreover, we may now state in 
addition that E C^([0, T], i7^(M")) and use Duhamel's formula again to show that also € 

ei,([0,r],iJi(M")) C W^^°^{[0,T],H\R")) as claimed. 

Step 2 (higher x- derivatives): We take the partial .x, -derivative on both sides in Equations (6) 
and (7) to obtain a similar differential equation for dx^u^ in the form 

n n 

dtdxjUe - i ^ dx^iCkedxkdxjUe) - tVedxjUe = i ^ dx^idxjCkedxkUe) + idx^VeUe + dxje =■ fie 

k=l k=l 

and the initial condition dxjUeiO) = dx^Qe- Thanks to (*) we have fu € L^([0, T], L^(R")) and 
dtfu G L^{[0,T],H'^{R")), which corresponds to conditions (7.5-6) of [11, Chapter XVIII, §7, 
Section 1, Theorem 1]. Thus we may apply (^i) and (8) with dx^Ue, dx^Qs, fis replacing u^, g^, 
fe, respectively, but with the same constants Ci and C2. Collecting the results for j = 1, . . . ,n 
we arrive at 

{S2) We e e([o,T],if2(M"))nei([o,T],L2(R")) 

and 

BiVseN: \\dx,dx,u,\\ <nVT sup ||«,(t)||^, = 0(£-^^) {e ^ 0, j,l = 1, . . . ,n). 

^ ^ ' " te[o,T] 

Similarly, taking higher partial x-derivatives in Equations (6) and (7) produces the same kind of 
differential equations and initial conditions for these higher derivatives, where we may always con- 
sider the already estimated lower order derivatives of as part of the right-hand side. Therefore 
we obtain successively for arbitrary m e N 

{Sm) u, e e([o,T],if'"(M"))ne^([o,T],if™-i(iR")) 

and 

Step 3 (t- derivatives): Wc first note that due to (S*,,,) we have ^{^^(O) G i/°°(K") and Equation 
(6) evaluated at t = provides moderate e-asymptotics of its Sobolev norms. Taking the partial 
t-derivative in Equation (6) now yields again the same kind of differential equation for dtU^: as 
above when interpreted in the following way: 

n n 

dtdtUe - i ^dx^iCkedx^dtUe) - iVedtUe = i ^dx^{dtCkedxkUe) + idtVeUe + dtfe =■ fu- 

fc=l fc=l 

By property (Sm) (now replacing (*) in a similar argument in Step 2) we deduce the required 
regularity conditions on fu and dtfis G L^([0, T], if~^(M")) to apply {Si) and (8) now with dtU^, 
dtUs{0), fu replacing u^, g^, fe, respectively. We obtain for a G Ng, \oi\ < 1, 
Ue e e\[0,T],H\R"))ne\[0,T],H-\M.")), 3Mi e N : ||ata>e|L2(M.x[o,T]) = 0{s-^^). 

Furthermore, inserting here the same procedure as in Step 1 now yields for any m € N and for all 

a G Nff, |a| < m, 

G e'{[0,T],H"\R"j) n e2([0,r],i/"-i(M")), 3M„ G N : = ©(e-^-). 

Thus in turn we may now deduce from the above differential equation for dtu^ that ^^^^(O) G 
H°°{M."^) holds with moderate spatial Sobolev norms. Hence we may now play the same game 
again with dfu,. etc. and finally arrive at the statements that for arbitrary d, m G N and for all 

a G Wq, \a\ < m, we have 

Ue G e'^([0,r],if™(K"))ne''+2([o,T],i7"-i(K")), 3M„<i e N : ||a^^5,"w,||^.(j,„^[o_^]) = Oie-^--). 
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Therefore we may conclude that (u^) is moderate net and its class Sz,2(M" x [0,T]) defines a 
solution u to the Cauchy problem (4-5). 

The proof of uniqueness requires to show negligibility of any solution (m^) assuming that the 
right-hand side (/e) and the initial data {g^) are negligible. But this follows successively from 
the corresponding variant of the energy estimate (8) applied at each step of the sequence of 
differentiated differential equations used in course of the existence proof. □ 



We note that in case of smooth coefficients an easy integration by parts argument shows that any 
solution to the Cauchy problem obtained from the variational mcrthod as in [11, Chapter XVIII, 
§7, Section 1]) also provides a solution in the sense of distributions. The following statement 
ensures in addition the coherence with the Colombeau generalized solution. 

Corollary 3.2. Let V and Ck (k = l,...,n) belong to C°°(f2T) n L°°{Qt) with bounded time 

derivatives of first- order, go G H^(M."), a,nd fo G C"^ ([0, T], iy^(M")). Let u denote the unique 
Colombeau generalized solution to the Cauchy problem (4-5), where g, f denote standard embed- 
dings of go, fo, respectively. Then up^ w, where w G C([0, T], is the unique distributional 

solution obtained from the variational method. 



Proof. Let (ug) e u be as in the proof of Theorem 3.1, then := Ue — w satisfies the following 
Cauchy problem 

n 

dtVe -^^dx^ (ckdxkVe) " iVVg = fe- .fo, We |t=0= 9e " 50- 

fe=l 

Moreover, we have the energy estimate (8), where now both Cu =: Ci and C2s ='■ C2 both are 
independent of e, in the form 

T 

sup Wveml^ < C2e^^(||5e - Qofm + [ (WMr) - fo{r)\\l. + ||9t/e(r) - dtfo{T)\\l-^) dr). 







Since the smoothing process via mollification ensures convergence to of the norms on the right- 
hand side of the above estimate we obtain even norm convergence of to w. □ 

Remark 3.3. To compare our main result with the evolution systems solution constructed in 
[12] for the case of Sobolev coefficients and spatial dimension n = 2 we drop the additional 
pseudodifferential aspects in that model. Then we obtain also coherence of the Colombeau concept 
with the functional analytic solution as we sketch in the following discussion: 

Let < r < 1 and y = 0, cin,C2o G ei([0, T], ff^'+HM^)), <7o e H'^i^^), fo G ei([0, T], ^^(R^)). 
Let V denote the unique iJ^-valuecl solution to the Cauchy problem corresponding to (4-5) with 
initial value go and right-hand side fo (and coefficients Cjo replacing Cj) according to [12, Theorem 
4.2]. Let 5, /, ci, C2 denote the embeddings of go, fo, cio, C20 into corresponding Colombeau 
spaces, respectively, and let u G Sl2(IR^ x [0,T]) be the unique Colombeau solution to (4-5). 
Subtracting the equations satisfied by a representative (u^) of u and v we deduce the following 
equations for the difference he := Ue — v 

2 2 

dthe - i ^ ^xfc {Ckedxkhs) = fe - fo + i^d^k ((Cfce " Cko)dxkV) =■ fe, h^ \t=0= 9e - 90, 

fe=l fe=l 
which imply an energy estimate of the form (8) for he and with suitably adapted initial value 

and right-hand side fe instead. As noted in the proof of Corollary 3.2 we have appropriate norm 
convergence of the standard regularizations ge, fe of the data go,fo- In addition, we now have to call 
on uniform boundedness of ||5tCe||j;,oo(][{2x[o t]) C^([0, T], iJ'"+^(M^))-norm convergence 

Cke - Cfco ^ (e ^ 0). Using the continuity of H'"+'^{R^) ■ H'^{M.^) H'"+^{R'^) we obtain A -> 
in ei([0,r],iJ'^+i(IR2)). Thus we finally obtain the convergence Ue ^ v in e{[0,T], H^(R^)), and 
in particular we have m ~ u. 
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Example 3.4. Consider the strictly positive square root of S represented by {^/Pe)ee]o,l]^ where 
p is a moUifier as in Proposition 2.1. (Note that wc obtain h,, = S * = in this case.) To 
simplify technical matters, assume in addition that p satisfies the conditions discussed in Remark 
2.2(ii) and moreover ^ G L^(IR"). (For example, any suitably normalized function of the form 

.T (1 + lajp)"™/^ with m > 2n would do.) Then we may consider the class g € Si,2(K"), given 
directly by {y/]h)ee]o,i] - as a square root of 5. 

We have g"^ « 5, but 5 « 0, which is easily seen by action on a test function ^p upon substituting 
y = x/e in / ^J Pe[x) (f{x)dx = e"/^ J y^p(y) tp{ey)dy and applying dominated convergence (thereby 
using that ^/p G L^). In essence, this effect has already been observed earlier in the generalized 
function model of ultrarelativistic Reissner-Nordstr0m fields in [33, Equations (15) and (17)]. 

The Cauchy problem (4-5) with generalized initial value g, right-hand side / = 0, constant coeffi- 
cients Cfe = 1 (fc = 1, . . . , n) and V = 0, written out for representatives then reads 

The solution is given by the action of the strongly continuous unitary group Ut ■— exp(itA) [t e M) 
of operators on Z/^(M") in the form Us{t,x) = {Ut^/p^){x). 

Let f e M and /x* denote the positive measure on M" with density function [^^(t, .)p for the 
Lebesgue measure. By unitarity of Ut we obtain 

4(K") = j \u,{t,x)\''dx = j{Ut^e){x)-{Ut^,){x)dx = J \^/^\''dx = J p,{x)dx = l, 

R" R" R" R" 

hence {/x* : t £ M, e g]0, 1]} is a family probability measures on R" and ||A*s||va,r = 1 (i € M, 
e e ]0, 1]) holds in the Banach space of finite measures. 

We claim that for any t the net {l-il)ee]o,i] converges to with respect to the vague topology 
on finite measures (cf. [3, §30]): Since ^/p^ € Z/^(IR") we obtain from the L^-L°°-estimate for 
the Schrodinger propagator ([31, §4.4, Theorem 1]) that [[^^(t, .)||^oo < ||\/p7|lz,i/(47r|i|)"/2 and 
therefore for any ip G Cc(IR.") 



|(4,V)I< I \u,{t,x)f\mx)\dx<{4n\t\)-^ II V-ILi 11x^711 



= (47r|tr"||^||^,||Vp||^,e"/2-^0 (e^O). 

The same obviously holds with -0 € ^(M"), hence also /x* — >■ in o5^'(R") (when t ^ 0), whereas 
(/i*)£g]o,i] can certainly not be weakly convergent as net of finite measures (i.e., in the weak-* 
topology in the dual of the space of bounded continuous functions on M") since (/x* , 1) = /U* (M") = 
1 7^ as e ^ (cf. also [3, Theorem 30.8]). 

Acknowledgement: The author wishes to express his sincere thanks to two anonymous referees. 
In particular, we thank the one for hints to additional references of related work, and we are 
deeply grateful to the second referee for patient comments and clarifications that lead to crucial 
corrections in Theorem 3.1 and its proof. 
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